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An analytical model of dispersive pulse propagation in semi-infinite cylinders due to transient 
axially symmetric end conditions has been experimentally investigated. Specifically, the dispersive 
propagation of the first axially symmetric longitudinal mode in thin wire waveguides, which have 
ends in butt contact with longitudinal piezoelectric ultrasonic transducers, is examined. The method 
allows for prediction of a propagated waveform given a measured source waveform, together with 
the material properties of the cylinder. Alternatively, the source waveform can be extracted from 
measurement of the propagated waveform. The material properties required for implementation of 
the pulse propagation model are determined using guided wave phase velocity measurements. Hard 
tempered aluminum 1100 and 304 stainless steel wires, with 127, 305, and 406 /im diam, were 
examined. In general, the drawn wires were found to behave as transversely isotropic media. 

PACS numbers: 43.35. Pt, 43.20.Jr 


INTRODUCTION 

Historically, pulse propagation in thin cylindrical 
waveguides has been of interest for applications including 
the material testing of wires and fibers, 1 3 ultrasonic delay 
lines, 4 medical ultrasonics, 5 and embedded sensor materials. h 
Propagation of ultrasonic waves in cylinders has been the 
subject of several investigations, both analytical and 
experimental. 7-20 In high-frequency and broadband applica- 
tions, these cylinders may exhibit considerable dispersion 
when an acoustic pulse propagates through the waveguide. 
Here, dispersion is the dependence of acoustic guided wave 
phase velocity on frequency, resulting in distortion of the 
waveform of a transient plane-wave pulse depending on the 
axial location. This pulse distortion can make the interpreta- 
tion of measured waveforms difficult. The present paper 
combines experimental measurements and analytical model- 
ing to yield a quantitative and integrated approach to address 
the problem of pulse dispersion in cylindrical waveguides. 
The method presented also has significant relevance to the 
problem of material characterization of these waveguides, 

I. PULSE PROPAGATION 

While exact 7 ' 9 and approximate 10,11 solutions have been 
used to study the axially symmetric guided wave phase ve- 
locity versus frequency relationship in solid cylindrical 
waveguides, and experimental dispersion measurements 
have corroborated the theory, 12-15 few studies have been 
done to assess the capability these solutions yield for mod- 
eling actual pulse propagation. The earliest analytical solu- 
tion, using exact theory, to the problem of axially symmetric 
pulse propagation in a cylindrical waveguide due to applica- 
tion of a transient end condition can be attributed to 
Skalak. 16 A special solution was obtained to the problem of a 
solid isotropic cylinder with a velocity imparted uniformly to 


the cylinder end in the form of a step function. Subsequently, 
Folk et al , 17 used exact theory to obtain a solution for a solid 
isotropic cylinder with arbitrary time harmonic mixed (com- 
bined displacement and stress) end boundary conditions. 
Moreover, an attempt was made by Folk et al . 1 7 to experi- 
mentally validate the solution in terms of its ability to model 
pulse propagation. Studies by Miklowitz and Nisewanger, 18 
Hseih and Kolsky, 19 and McNiven and Mengi 20 have at- 
tempted to experimentally validate propagation models 
based on approximate theories. Only qualitative agreement 
between theoretical and measured waveforms was obtained 
in these studies. Furthermore, the pulses examined in all of 
these studies were dominated by lower-frequency propaga- 
tion of the first longitudinal mode. The pulses, resulting from 
either ideal or modified (some finite rise time) step loading, 
essentially consist of a low-frequency leading head and a 
higher-frequency dispersive tail. While the current study also 
addresses first longitudinal mode propagation, the transient 
end conditions applied give rise to more varied dispersive 
pulse distortion. The frequency range considered extends 
from low frequency, where the fundamental tends toward 
nondispersiveness, to approaching the second longitudinal 
mode cutoff, where the fundamental becomes highly disper- 
sive. An integrated analytical and experimental study of iso- 
tropic hollow cylinders using the same numerical method 
and similar measurement techniques 21 involves the propaga- 
tion of higher-order longitudinal modes. 


II. PROPERTY MEASUREMENT 

Because the dispersion of the guided waves depends on 
the elastic properties, cylindrical guided wave propagation 
has been used for the material characterization and testing of 
wires and fibers. Meeker and Meitzler 22 suggested several 
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methods for determining material properties by measuring 
acoustic guided wave velocities, and a thorough review of 
the subject has been given by Curtis. 1 

Principally, velocity measurements of the first torsional 
mode and the long wavelength limit of the first longitudinal 
mode have been used to obtain the shear and Young’s 
moduli, respectively, for isotropic cylinders. 1 2 The first axi- 
symmetric torsional mode is nondispersive, propagating at 
the shear wave velocity (G/p) [ 2 . In the long wavelength 
(low-frequency) limit, the velocity of the first longitudinal 
axisymmetric mode tends toward the bar velocity (E/p) 12 . 
As the limit is approached, group and phase velocities tend 
toward a constant and equal value, indicating nondispersive 
behavior. 

However, the drawing process in the manufacture of 
wires and fibers can affect the material microstructure, fre- 
quently causing the material to be anisotropic. 1 Transverse 
isotropy, with five independent elastic constants, is consid- 
ered to better represent the state of the material than isotropy. 
The five independent elastic properties for transversely iso- 
tropic media are taken as £ LL , E jt, G LI , (7^, v LT , where 
E, G\ and v are Young's moduli, shear moduli, and poisson 
ratio, respectively, L is longitudinal, and T is transverse. 23 
Moreover, the propagation of the fundamental longitudinal 
mode, within the frequency range considered, does not de- 
pend on the transversely isotropic property to which the low- 
est torsional mode is sensitive, G ir . Rather, the strongest 
dependence is on E LL and y LT . 

Longitudinal mode measurement to obtain E LL can be 
problematic. When the velocity of the long wavelength limit 
is used, a transient pulse with sufficient low-frequency con- 
tent is propagated and its first arrival measured. Since the 
longest wavelengths of the first longitudinal mode propagate 
with greatest velocity, this time-of-flight measurement yields 
the velocity of the long wavelength limit. However, the lead- 
ing head of a dispersive pulse often rises very gradually mak- 
ing accurate assessment of when the first arrival occurs dif- 
ficult, as discussed by Curtis, 1 and evident in the low- 
frequency experimental waveforms of other studies. 17-20 
Consequently, the measurement may be subject to some er- 
ror. Instead, the time shift between the peak of the wave 
envelope of a propagated broadband pulse is often matched 
with that of the initial pulse. 2 However, the frequency and 
wavelength of such a velocity measurement is not at all 
clear. 24 

Alternatively, phase or group velocity of particular cy- 
lindrical guided wave modes can be measured over a range 
of frequency. Material properties are then varied to obtain 
correspondence between theoretical velocities and the disper- 
sion measurements. As noted by Meeker and Meitzler, 2 " dis- 
persion measurements of the first axially symmetric longitu- 
dinal mode are sufficient to find both elastic constants for 
isotropic media. As will be illustrated here, combination of 
such measurement with measures of first torsional mode ve- 
locity yields E LL , v u , and G lT for transversely isotropic 
media. The extent of the material anisotropy can then be 
gauged. 

Phase velocity measurement methods are preferred to 
group velocity measurement methods since the wave propa- 
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FIG. 1. Experimental apparatus. 

gation is generally dispersive. Inherent difficulty in group 
velocity measurements for dispersive media has been 
noted. 24,25 The problem stems from the distortion of the 
propagated waveform relative to the initially generated 
waveform and the resulting ambiguity in correlating the two 
different waveforms, and also, in determining the frequency 
for which the group velocity is being measured. Correspond- 
ingly, precision and accuracy for phase velocity measure- 
ment can be expected to exceed that obtained for group ve- 
locity measurement for dispersive pulse propagation. The 
scatter of group velocity measurements which commonly oc- 
curs for propagation in dispersive media can be seen in a 
recent study of sound in wires by Nicholson etui 15 The 
precision of the measurement and smoothness of the curve 
determined by these measurements is critical to accurate and 
repeatable material property characterization. The elastic 
properties determined from the phase velocity measurements 
are used here as input to the pulse propagation model. 

III. WAVE EXCITATION AND MEASUREMENT 

The method by which elastic pulses are launched into 
the wire specimens and measured after propagating the 
length of the sample is illustrated in Fig. 1. Similar tech- 
niques have been used in other studies of ultrasound propa- 
gation in wires. 15 As shown in Fig. 1, matched pairs of pi- 
ezoelectric longitudinal or shear wave transducers are placed 
in butt contact with each end of the wire and operate in pitch 
catch mode. For this study, 5 MHz center frequency pana- 
metrics VI 09 (longitudinal) and VI 56 (shear) transducers 
were used. Glycerin was used as couplant to aid in sound 
transmission across the transducer/wire interface. The source 
or excitation transducer was provided a gated sine wave or 
toneburst electrical signal. Because of the thinness of the 
wires, the efficiency at which sound could be introduced into 
the waveguides was very low. Consequently, gain in the re- 
ceiving system between 60 and 100 dB was required and 
signal averaging between 300 and 600 data records was per- 
formed to improve signal to noise. Signals were captured and 
averaged on a digital oscilloscope and processed on a PC. 

With longitudinal wave transducers in use, the pulse 
launched in the medium would be composed primarily of the 
first axially symmetric longitudinal mode. However, depend- 
ing on factors such as the wire thickness, excitation fre- 

Madaras et al.: Dispersive pulse propagation in wire 253 



Measurement and modeling of dispersive pulse propagation 
in drawn wire waveguides 

Eric I. Madaras 

NASA Langley Research Center , ; Hampton, Virginia 236X1-0001 

Thomas W. Kohl and Wayne P. Rogers 

University of Colorado, Department of Mechanical Engineering, Boulder, Colorado 80309-0427 
(Received 20 October 1993; accepted for publication 12 July 1994) 

An analytical model of dispersive pulse propagation in semi-infinite cylinders due to transient 
axially symmetric end conditions has been experimentally investigated. Specifically, the dispersive 
propagation of the first axially symmetric longitudinal mode in thin wire waveguides, which have 
ends in butt contact with longitudinal piezoelectric ultrasonic transducers, is examined. The method 
allows for prediction of a propagated waveform given a measured source waveform, together with 
the material properties of the cylinder. Alternatively, the source waveform can be extracted from 
measurement of the propagated waveform. The material properties required for implementation of 
the pulse propagation model are determined using guided wave phase velocity measurements. Hard 
tempered aluminum 1100 and 304 stainless steel wires, with 127, 305, and 406 /x m diam, were 
examined. In general, the drawn wires were found to behave as transversely isotropic media. 

PACS numbers: 43.35.Pt, 43.20Jr 


INTRODUCTION 

Historically, pulse propagation in thin cylindrical 
waveguides has been of interest for applications including 
the material testing of wires and fibers, 1 ' 3 ultrasonic delay 
lines, 4 medical ultrasonics, 5 and embedded sensor materials. 6 
Propagation of ultrasonic waves in cylinders has been the 
subject of several investigations, both analytical and 
experimental. 7 ' 20 In high-frequency and broadband applica- 
tions, these cylinders may exhibit considerable dispersion 
when an acoustic pulse propagates through the waveguide. 
Here, dispersion is the dependence of acoustic guided wave 
phase velocity on frequency, resulting in distortion of the 
waveform of a transient plane-wave pulse depending on the 
axial location. This pulse distortion can make the interpreta- 
tion of measured waveforms difficult. The present paper 
combines experimental measurements and analytical model- 
ing to yield a quantitative and integrated approach to address 
the problem of pulse dispersion in cylindrical waveguides. 
The method presented also has significant relevance to the 
problem of material characterization of these waveguides. 

I. PULSE PROPAGATION 

While exact 7 ' 9 and approximate 10,11 solutions have been 
used to study the axially symmetric guided wave phase ve- 
locity versus frequency relationship in solid cylindrical 
waveguides, and experimental dispersion measurements 
have corroborated the theory, 12 ' 15 few studies have been 
done to assess the capability these solutions yield for mod- 
eling actual pulse propagation. The earliest analytical solu- 
tion, using exact theory, to the problem of axially symmetric 
pulse propagation in a cylindrical waveguide due to applica- 
tion of a transient end condition can be attributed to 
Skalak. 16 A special solution was obtained to the problem of a 
solid isotropic cylinder with a velocity imparted uniformly to 


the cylinder end in the form of a step function. Subsequently, 
Folk et al} 1 used exact theory to obtain a solution for a solid 
isotropic cylinder with arbitrary time harmonic mixed (com- 
bined displacement and stress) end boundary conditions. 
Moreover, an attempt was made by Folk et al} 1 to experi- 
mentally validate the solution in terms of its ability to model 
pulse propagation. Studies by Miklowitz and Nisewanger, 18 
Hseih and Kolsky, 19 and McNiven and Mengi 20 have at- 
tempted to experimentally validate propagation models 
based on approximate theories. Only qualitative agreement 
between theoretical and measured waveforms was obtained 
in these studies. Furthermore, the pulses examined in all of 
these studies were dominated by lower-frequency propaga- 
tion of the first longitudinal mode. The pulses, resulting from 
either ideal or modified (some finite rise time) step loading, 
essentially consist of a low-frequency leading head and a 
higher-frequency dispersive tail. While the current study also 
addresses first longitudinal mode propagation, the transient 
end conditions applied give rise to more varied dispersive 
pulse distortion. The frequency range considered extends 
from low frequency, where the fundamental tends toward 
nondispersiveness, to approaching the second longitudinal 
mode cutoff, where the fundamental becomes highly disper- 
sive. An integrated analytical and experimental study of iso- 
tropic hollow cylinders using the same numerical method 
and similar measurement techniques 21 involves the propaga- 
tion of higher-order longitudinal modes. 


II. PROPERTY MEASUREMENT 

Because the dispersion of the guided waves depends on 
the elastic properties, cylindrical guided wave propagation 
has been used for the material characterization and testing of 
wires and fibers. Meeker and Meitzler 22 suggested several 
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methods for determining material properties by measuring 
acoustic guided wave velocities, and a thorough review of 
the subject has been given by Curtis. 1 

Principally, velocity measurements of the first torsional 
mode and the long wavelength limit of the first longitudinal 
mode have been used to obtain the shear and Young's 
moduli, respectively, for isotropic cylinders. 1,2 The first axi- 
symmetric torsional mode is nondispersive, propagating at 
the shear wave velocity (G/p) 1 2 . In the long wavelength 
(low-frequency) limit, the velocity of the first longitudinal 
axisymmetric mode tends toward the bar velocity (£7p) I/2 . 
As the limit is approached, group and phase velocities tend 
toward a constant and equal value, indicating nondispersive 
behavior. 

However, the drawing process in the manufacture of 
wires and fibers can affect the material microstructure, fre- 
quently causing the material to be anisotropic. 1 Transverse 
isotropy, with five independent elastic constants, is consid- 
ered to better represent the state of the material than isotropy. 
The five independent elastic properties for transversely iso- 
tropic media are taken as £ LL , E rr , ^rr» where 
E , G, and v are Young’s moduli, shear moduli, and poisson 
ratio, respectively, L is longitudinal, and T is transverse. 23 
Moreover, the propagation of the fundamental longitudinal 
mode, within the frequency range considered, does not de- 
pend on the transversely isotropic property to which the low- 
est torsional mode is sensitive, G L t- Rather, the strongest 
dependence is on E LL and y LT . 

Longitudinal mode measurement to obtain E LL can be 
problematic. When the velocity of the long wavelength limit 
is used, a transient pulse with sufficient low-frequency con- 
tent is propagated and its first arrival measured. Since the 
longest wavelengths of the first longitudinal mode propagate 
with greatest velocity, this time-of-flight measurement yields 
the velocity of the long wavelength limit. However, the lead- 
ing head of a dispersive pulse often rises very gradually mak- 
ing accurate assessment of when the first arrival occurs dif- 
ficult, as discussed by Curtis, 1 and evident in the low- 
frequency experimental waveforms of other studies. 17-20 
Consequently, the measurement may be subject to some er- 
ror. Instead, the time shift between the peak of the wave 
envelope of a propagated broadband pulse is often matched 
with that of the initial pulse. 2 However, the frequency and 
wavelength of such a velocity measurement is not at all 
clear. 24 

Alternatively, phase or group velocity of particular cy- 
lindrical guided wave modes can be measured over a range 
of frequency. Material properties are then varied to obtain 
correspondence between theoretical velocities and the disper- 
sion measurements. As noted by Meeker and Meitzler,^ 2 dis- 
persion measurements of the first axially symmetric longitu- 
dinal mode are sufficient to find both elastic constants for 
isotropic media. As will be illustrated here, combination of 
such measurement with measures of first torsional mode ve- 
locity yields E LL , v u , and G u for transversely isotropic 
media. The extent of the material anisotropy can then be 
gauged. 

Phase velocity measurement methods are preferred to 
group velocity measurement methods since the wave propa- 
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FIG. 1. Experimental apparatus. 

gation is generally dispersive. Inherent difficulty in group 
velocity measurements for dispersive media has been 
noted. 24,25 The problem stems from the distortion of the 
propagated waveform relative to the initially generated 
waveform and the resulting ambiguity in correlating the two 
different waveforms, and also, in determining the frequency 
for which the group velocity is being measured. Correspond- 
ingly, precision and accuracy for phase velocity measure- 
ment can be expected to exceed that obtained for group ve- 
locity measurement for dispersive pulse propagation. The 
scatter of group velocity measurements which commonly oc- 
curs for propagation in dispersive media can be seen in a 
recent study of sound in wires by Nicholson et al . 15 The 
precision of the measurement and smoothness of the curve 
determined by these measurements is critical to accurate and 
repeatable material property characterization. The elastic 
properties determined from the phase velocity measurements 
are used here as input to the pulse propagation model. 

III. WAVE EXCITATION AND MEASUREMENT 

The method by which elastic pulses are launched into 
the wire specimens and measured after propagating the 
length of the sample is illustrated in Fig. 1. Similar tech- 
niques have been used in other studies of ultrasound propa- 
gation in wires. 15 As shown in Fig. 1, matched pairs of pi- 
ezoelectric longitudinal or shear wave transducers are placed 
in butt contact with each end of the wire and operate in pitch 
catch mode. For this study, 5 MHz center frequency pana- 
metrics V109 (longitudinal) and V156 (shear) transducers 
were used. Glycerin was used as couplant to aid in sound 
transmission across the transducer/wire interface. The source 
or excitation transducer was provided a gated sine wave or 
toneburst electrical signal. Because of the thinness of the 
wires, the efficiency at which sound could be introduced into 
the waveguides was very low. Consequently, gain in the re- 
ceiving system between 60 and 100 dB was required and 
signal averaging between 300 and 600 data records was per- 
formed to improve signal to noise. Signals were captured and 
averaged on a digital oscilloscope and processed on a PC. 

With longitudinal wave transducers in use, the pulse 
launched in the medium would be composed primarily of the 
first axially symmetric longitudinal mode. However, depend- 
ing on factors such as the wire thickness, excitation fre- 
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FIG. 2. Measured propagated pulse in a 24.76-cm-long, 12-mil-diam alumi 
num 1100 wire due to excitation of longitudinal piezoelectric source trans- 
ducer with 1-MHz (top) and 5-MHz (bottom) tonebursts. 


FIG. 3. Measured propagated pulse in a 25.22-cm-long, 5-mil-diam 304 
stainless steel wire due of excitation to shear piezoelectric source transducer 
with 700-kHz (top) and 1. 5-MHz (bottom) tonebursts. 


quency, couplant effects, and obviously the contact between 
source transducer and wire, the first-order flexural mode was 
also excited. In general, the flexural mode propagates much 
more slowly than the axisymmetric longitudinal mode. Con- 
sequently, the waveform due to each of the two modes can be 
temporally separated after some length of propagation. How- 
ever, with increasing frequency, the dispersion of the wave- 
guide is such that the first longitudinal mode is observed to 
propagate more slowly while the flexural mode is observed 
to propagate faster. Propagation of higher-frequency harmon- 
ics may result in overlap and interference between axially 
symmetric longitudinal and nonsymmetric flexural waveform 
components. The longitudinal and flexural components of 
typical measured waveforms for bursts with increasing cen- 
ter frequency are shown in Fig. 2. Longitudinal wave trans- 
ducers are used for generating and receiving the wave. The 
waveform component associated with the flexural mode was 
eliminated through windowing for both the phase velocity 
measurements and the extraction of the source signal from a 
measured received waveform. 

Torsional modes were excited by using shear wave 
transducers and intentionally placing only an edge of the 
wire end into contact with the source transducers face. As a 
result of such contact, both longitudinal and torsional axi- 


symmetric modes and at least first-order flexural modes were 
excited. By varying the frequency of the toneburst excitation 
of the shear wave transducer, and observing which wave 
packet remained fixed with regard to arrival and did not dis- 
tort relative to the source waveform (nondispersive), the tor- 
sional wave packet could be identified. Torsional wave pack- 
ets are shown for bursts of 0.7 and 1.5 MHz in Fig. 3. Since 
dispersive first-order flexural waves are also excited, caution 
must be used as flexural waves will often propagate at simi- 
lar velocities relative to torsional waves. Consequently there 
is overlap of the waveforms over broad ranges of frequency. 
In Fig. 3, the faster propagation of the flexural mode when a 
higher-frequency burst was used, results in the beginning of 
interference between the nondispersive torsional mode and 
the highly dispersive flexural mode. 

Source measurements were performed by placing the 
matched transducers into face to face contact without the 
intervening wire. It should be noted that source measurement 
by this method is imperfect since it neglects all transducer/ 
wire interface effects. In particular, some variation in the 
behavior of the source transducer when electrically excited in 
unloaded (with wire contacting transducer face) as opposed 
to loaded (matched transducer face-to-face contact) condition 
could be anticipated. The effect of multiple reflections be- 
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TABLE I. List of measured material properties. 


Wire type (diameter 
and density) 

£ L l (GPa) 

^LT 

Glt (GPa) 

G lso (GPa) 

304 stainless steel 





406 /xm (16 mil), 
7910 kg/m 3 

211 

0.30 

74.5 

81 

305 \x m (12 mil), 
7852 kg/m 3 
(annealed) 

208 

0.30 

76.0 

80 

127 ptm (5 mil), 
7913 kg/m 3 

225 

0.35 

71.1 

82 

aluminum 1100 





406 /xm (16 mil), 
2679 kg/m 3 

72.8 

0.36 

24.4 

27 

305 /xm (12 mil), 
2668 kg/m 3 

68.6 

0.36 

25.3 

25 

127 /xm (5 mil), 
2742 kg/m 3 

71.0 

0.40 

26.0 

25 


increases. Most real dispersive pulses will have regions of 
frequency exhibiting nonphysically realizable phase (<)>) in- 
formation (i.e., corresponding to negative group velocity), 
often concomitant with low-magnitude complex amplitude 
A(w). With a high density of phase jumps and even small 
regions of inaccurate phase information, the phase wraps be- 
come impossible to accurately track. 

As noted previously, curve fitting of measured phase ve- 
locities of the first longitudinal axisymmetric mode through 
variation of material properties allows complete material 
property characterization for waveguides of isotropic 
medial 2 The same method may be applied to transversely 
isotropic cylinders provided that E ll an d J'lt are token as the 
independent properties to be varied in obtaining the curve fit. 
Phase velocity measurements from low frequency to fre- 
quencies evidencing significant dispersion (approaching the 
cutoff of the second longitudinal mode) are required for ac- 
curate determination of both material properties. 

The lowest torsional axisymmetric mode is nondisper- 
sive, for both transversely isotropic and isotropic cylinders. 
This mode is sensitive only to G L j for transverse isotropy. 
Once the torsional wavepacket is identified, its velocity is 
measured by cross-correlation." 8 

Table I summarizes the properties obtained in this man- 
ner for six different wire types investigated. The final column 
for each wire gives the shear modulus which would be re- 
quired to support the conclusion of isotropy given the values 
of E ll and v LT determined by longitudinal mode measure- 
ment. Comparing this value of shear modulus with that ob- 
tained by torsional mode measurement is a measure of the 
material anisotropy and supports the conclusion that trans- 
verse isotropy better corresponds to the material behavior 
than isotropy, particularly for the steel wires. However, there 
is good agreement between solutions for the dispersion rela- 
tions using either Pochhammer’s exact solution for isotropic 
media taking E and v as independent properties, or using the 


numerical method described above and the complete set of 
five independent properties. The implication is that determi- 
nation of G IT was only of interest in verifying the anisotropy 
of the media, but this anisotropy has no significant effect on 
the first longitudinal mode behavior. Since the first longitu- 
dinal mode propagation was not sensitive to the exact values 
for Ejj and G^, these were set equal to E LL and G LT , 
respectively. 

Figures 4 and 5 show measured longitudinal mode ve- 
locities and the curve fits obtained, for 16-mil (406 yitm) and 
5-mil (127 yam) diameter stainless steel wires, respectively. 
The entire frequency range for which data was taken was not 
necessarily used for each wire. The solid circles indicate 
phase velocity measurements at discrete frequencies used in 
the curve fit shown as a solid line. In Fig. 5, additional mea- 
surements not used in the curve fit are shown as squares. 



FIG. 4. Curve fit (solid line) of measured phase velocity (solid circles) to 
determine elastic properties for a 16-mil-diam 304 stainless steel wire. Ana- 
lytical curves due to 0.5% increase in £ LL (dotted line) and 5% increase in 
v (dot-dash line) are also shown. 
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Frequency (MHz) 



FIG. 5. Curve fit (solid line) of measured phase velocity (solid circles) to 
determine clastic properties for a 5-mil-diam 304 stainless steel wire. Ana- 
lytical curves due to 0.5% increase in E LI (dotted line) and 5% increase in 
v (dot-dash line) are also shown. 


Additional information regarding error analysis and accuracy 
is also conveyed and requires some explanation. 

The appearance of the measured dispersion at lower fre- 
quencies in Fig. 5 does not correspond to the first mode 
dispersion for any elastic cylinder. Such measurement 
anomaly was much more significant for the smaller diameter 
(127 yum) wires tested. The frequency dependences of trans- 
mission through both of the transducer/wire interfaces, 
which has been neglected, may plausibly be influenced by 
the wire diameter. 

In general, several other factors may contribute to mea- 
surement error. First, the source waveform may be inad- 
equately characterized, due to the difference in source trans- 
ducer loadings between source measurement and 
transmission into the wire specimen. Furthermore, at the re- 
ceiving end of the wire, some of the wave incident on the 
transducer is reflected, which may also give rise to some 
phase shift. The frequency difference between lower- 
measurement frequencies and the resonant frequency of the 
transducer (5 MHz) may exacerbate these problems. 

If any resulting phase shift error at the pulse center fre- 
quency does not exceed some maximum value, say it/2, and 
the determination of M remains unambiguous, it is clear from 
Eq. (6a) that the phase velocity measurement becomes more 
accurate as frequency increases. The actual phase (<f>) mea- 
surement becomes increasingly small relative to the quantity 
27 tN in the denominator of Eq. (6a). The error bars of Figs. 
4 and 5 indicate the change in measured phase velocity 
which would result from a ±tt/2 shift in <f>. 

Because of the apparent systematic error which occurred 
at lower frequencies for the 5-mil steel wire (Fig. 5), not all 
of the dispersive data available was used in the curve fit. 
That the analytical curve fit falls within the range of the error 
bars of the unused measurements for a 7t/2 maximum error 
in <p, indicates that the phase shift error corresponds to less 
than a quarter of a wavelength. For the 16-mil diameter 
stainless steel wire (Fig. 4), measurements beyond 5 MHz 
were not used because of overlap and interference between 
longitudinal and flexural waves, as discussed previously. 

The dashed lines of Figs. 4 and 5 show the alteration of 


o 

> 


CO 

c 

o> 

oo 

<D 

> 

’55 

o 

a> 

QC 


1.5 
1.0 
0.5 
0.0 
- 0.5 
- 1.0 
- 1.5 

47.0 51.0 55.0 59.0 63.0 


/ 1 / 

> \ I 

(1 (1 

j* 

1 

1 1 — 

l!'- 

i . . i i 


(b) 


Time (microseconds) 


FIG. 6. Pulse propagation of a 1-MHz toneburst in 25.27-cm-long, 16-mil- 
diam 304 stainless steel wire; (a) measured (solid line) and analytically 
extracted (dashed line) source waveforms, (b) measured (solid line) and 
analytically predicted (dashed line) receiver waveforms. 


the analytically determined curves due to 0.5% and 5% 
variations in E LL and y LT , respectively. In comparing Figs. 4 
and 5, specifically the dashed and solid curves in each figure, 
it is apparent that in order to determine v LJ accurately, data 
associated with the increased dispersion at higher frequency, 
approaching the second longitudinal mode cutoff, must be 
available. Consequently, accuracy of the determination of 
v LT for the thicker wire was judged to be improved (i.e., the 
second mode cutoff is at much higher frequency for the thin- 
ner wire). The dispersion of the fundamental mode, as fre- 
quency and wire diameter to wavelength ratio increase, has 
been associated with the significance of radial inertia. 10,11 A 
conservative accuracy estimate of better than 0.5% for E LL 
and 5% for v LT is claimed when adequate dispersion mea- 
surements are available. 

VI. COMPARISON OF EXPERIMENTALLY MEASURED 
AND MODELED PULSE PROPAGATION 

Comparison of experimentally measured receiver wave- 
forms with analytically predicted responses, due to various 
measured sources, for the 16-mil steel wire, are shown in 
Figs. 6(b), 7(b), 8(b), and 9(b). The measured sources engen- 
dering these responses are shown in Figs. 6(a), 7(a), 8(a), and 
9(a) together with the source waveform extracted from the 
aforementioned measured receiver waveforms. The analyti- 
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